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Abstract 

We present a new mechanism for transferring a pre-existing lepton or baryon asymmetry to a dark matter 
asymmetry that relies on mass mixing which is dynamically induced in the early universe. Such mixing 
can succeed with only generic scales and operators and can give rise to distinctive relationships between 
the asymmetries in the two sectors. The mixing eliminates the need for the type of additional higher- 
dimensional operators that are inherent to many current asymmetric dark matter models. We consider 
several implementations of this idea. In one model, mass mixing is temporarily induced during a two-stage 
electroweak phase transition in a two Higgs doublet model. In the other class of models, mass mixing is 
induced by large field vacuum expectation values at high temperatures - either moduli fields or even more 
generic kinetic terms. Mass mixing models of this type can readily accommodate asymmetric dark matter 
masses ranging from 1 GeV to 100 TeV and expand the scope of possible relationships between the dark 
and visible sectors in such models. 
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1 Introduction 

Asymmetric dark matter models p] have the potential to explain the fact that the energy in dark matter and 
ordinary matter are notably comparable. Suggested explanations involve the co-generation of an asymmetry 
in both dark matter and baryonic sectors by the decay of heavy fields [2J [3] or the Affleck-Dine mechanism 
[4]. Many others involve higher-dimensional operators that transfer asymmetries between the two, as was first 
explored in [5] and more recently in p] |6j [3 HJ . Although it seems reasonable that dark matter number and 
baryon or lcpton number might not be exactly preserved in the early universe, most such models require new 
ingredients (scales or fields) to accommodate higher-dimensional operators that lead to requisite relationships 
between the dark matter and ordinary matter densities. Although conceptually compelling, these models leave 
open the question of whether we really expect the operators that violate baryon or lepton and dark matter 
numbers to exactly the right degree to really exist. Much more convincing would be models where there is 
mass mixing in the early universe arising from renormalizable interactions or operators suppressed only by 
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the Planck scale. Ultimately, it is experiments that will hopefully shed light on the precise nature of the dark 
sector and its connection to the Standard Model. 

In this paper, we show that it is possible that lepton (and hence baryon) number and dark matter number 
are not independent in the early universe. Mass mixing between dark-matter-number-carrying and lepton- 
number-carrying particles allows asymmetries to be established for both dark matter and lepton numbers. This 
mixing, of course, must turn off at later times when we know that the two symmetries must be independently 
preserved. The mechanism therefore requires fields whose value at high temperature differ from their low- 
temperature values so that in the late universe, dark matter and lepton numbers are independent symmetries 
and separately conserved. 

The mass parameters in the models we present are linked to scales we know should be associated with new 
fields or interactions, namely the Planck scale and the electroweak scale (although in some cases we invoke also 
leptogenesis or grand unification scales). Furthermore, mass mixing can allow for novel relationships between 
the asymmetries in the dark and visible sectors when asymmetry transfer processes are out of equilibrium or 
the mass mixing shuts off rapidly in a first-order phase transition, allowing for deviations from the often-cited 
"predicted" value of the dark matter mass of ~ 5 GeV. 

In our models, the mass mixing is active during a period of time after a fundamental B — L asymmetry 
is generated but vanishes at late times to ensure stability of the dark matter. This is the opposite of the 
usual chronology of phase transitions, where finite-temperature corrections keep the system in the symmetric 
phase at early times and a vacuum expectation value (VEV) arises at later times. There are at least two 
situations when symmetries are broken in the early universe and restored at late times. One involves multiple 
fields where one field initially stabilizes another field away from its minimum but no longer does so below a 
critical value [5]. Two stage phase transitions in a two Higgs model were first discussed in [TO]. A similar 
idea is exploited in hybrid inflation models, jllj . but the fields in our models are confined to the minima of 
the potential and involve no slow roll. We present models with multiple-stage phase transitions in Section [2] 
The second possibility is temperature-dependent background fields in the early universe that turn off as the 
universe cools. Examples include moduli field VEVs and the background energy density in the early universe. 
When these background fields are active, Planck-suppressed higher-dimensional operators can induce mass 
mixing at early times. We present both supcrsymmetric and non-supersymmetric models in Section [3] 

2 Two-stage phase transition 
2.1 Overview 

Canonical examples of phase transitions in the early universe involve a transition from the symmetric phase 
at high temperatures to a broken phase at low temperatures. In Weinberg's pioneering work on symmetry 
breaking at finite temperature [9], examples were presented where the reverse is true: a symmetry broken at 
high temperatures is restored at later times. Typically, there exist multiple scalar fields and multiple stages 
to such phase transitions. In this paper, we first discuss one simple example of a two Higgs doublet model 
originally considered in the context of electroweak baryogenesis pjj], although our results can be generalized 
in a relatively straightforward manner. The models in the following section are more generic, but have slightly 
more technically complicated dynamics and work for a more restricted parameter range. 

In the two Higgs set-up, a B — L asymmetry is established at temperatures above a few hundred GeV. The 
precise mechanism is irrelevant to our analysis. At high temperatures T 3> 100 GeV, the universe is in an 
electroweak symmetric phase. As the universe cools to T ~ 100 GeV, the electroweak phase transition occurs, 
proceeding in two stages. In the first stage, a non-Standard Model scalar $ acquires a VEV while the Higgs 
remains confined to the origin. When $ gets a VEV, it induces mass mixing between the dark sector and the 
left-handed Standard Model leptons. During this period, the lepton asymmetry gets shifted into an asymmetry 
involving the high-temperature mass eigenstates, which are linear combinations of the zero-temperature dark 
matter and lepton fields. 

At later time, a minimum develops along the Standard Model Higgs H direction with a lower energy than 
the initial vacuum. In the new vacuum, the Standard Model Higgs develops a VEV but $ is confined to the 
origin by the potential induced by the large Standard Model Higgs VEV. Bubbles of the new vacuum form in 
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a first-order phase transition. Because this phase transition occurs very rapidly, mass mixing shuts off before 
thermal interactions can adjust the asymmetry. This allows for non-trivial mixing of the asymmetry between 
the dark and visible sectors. In one scenario, the mixing can populate heavy dark sector fields that have 
been thermally suppressed in the previous stage, in which case the relationship between the final dark matter 
and lepton asymmetries is determined by the mixing angle at the time of the first-order phase transition. In 
another scenario, it is determined purely by thermal factors. 

In this two Higgs model, the second Higgs plays two important roles. First, it allows a field value that turns 
on and off before the fields enter the true electroweak-symmetry-breaking vacuum with the zero temperature 
Higgs VEV. Second, the <J> scalar is initially trapped far from the origin and allows for a first-order phase 
transition so the asymmetry in the dark matter gets established immediately. 

We propose models with 0(1) couplings and mass scales in the theory that are connected to the weak 
scale. We do not attempt to solve the hierarchy problem, but presume the model can be embedded in models 
that do have solutions (for instance, in supersymmetric models in which both scalars get their masses from 
supersymmetry breaking) . 



2.2 Asymmetry transfer during phase transition 

The mechanism of asymmetry transfer by mass mixing is largely independent of the specific details of the two 
Higgs model, and we postpone a discussion of the field content, interactions, and phase transition to Section 
12.31 and appendices [3] and |B] For now, we consider general models with a doublet scalar $ with a Yukawa 
coupling yx to a vector-like dark matter particle, A, and a left-handed lepton L. During the intermediate 
stage of the electroweak phase transition, the real part of the scalar i> has a VEV that induces mass mixing 
between A and L, 

m X L{T) = y x v$(T). (1) 

If $ does not acquire a VEV, an unbroken U(l)x-* symmetry would persist through the entire universe's 
evolution. At late times, after $ had decayed away entirely, the solution to the chemical potential relations 
would be /i$ = fix = 0, and there would be no late time X asymmetry. Thus, the mass mixing stage is vital 
to ensuring a non-zero asymmetry at late times. 

Although the model contains three generations of X and L, for simplicity and to retain only the essential 
elements of the mass mixing transfer mechanism, we consider only the case where one flavor of L and X have 
a large mixing. This could be obtained, for example, by hierarchical Yukawa couplings of the same form as in 
the Standard Model. The discussion can be easily generalized to the case of arbitrary mixings between flavor 
states. 

The vector-like dark matter has a Dirac mass mx- The mass matrix in the (X, X, L°) basis, including the 
dynamically induced mixing term mxL{T) is 

(0 m x mxh{T) 
m x 
m XL (T) 



which gives two degenerate states with mass M = m 2 x + m XL and a massless state. The mass basis states 
are 

\L') = -sin0LY) +cos(9|L°), (3) 
| A') = cos(9|A) + sin0|L°), (4) 

(5) 

with the \L') state being the massless state and the |A), |A) states having mass M. The mixing angle is 

• a mxL ia\ 

V m X + m XL 
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In the interaction basis, the rate of asymmetry transfer between L and A is approximately |12j 



/ A/2 \ 

r L ^=r sin 2 2#sin 2 ^— J, (7) 

which is explained in Section r3.2l This formula says that the net transfer rate is the oscillation rate that occurs 
before a thermal interaction, multiplied by the thermal interaction rate. Tq is the rate at which scattering 
events occur, averaged between A and L. The system is in chemical equilibrium for the mass scales and 
couplings considered in this section if y x > 10~ 6 . 

The gauge interaction of L° with charged leptons appears in the mass basis as 

/W - l(L-) f r(l-7 5 )W- (coseL' + sinex') . (8) 

The rate of interactions of X' is sin 2 9 times the usual weak interaction rate. Since the Hubble constant is 
very small at the time of the electroweak phase transition, for any non-negligible value of 9, X' is in thermal 
equilibrium due to gauge interactions. Since L~ carries lepton number of +1, X' and L' carry lepton number 
of +1 as well. The Yukawa interaction 

C D y x $ XL = y x $ X(- sin 9X' + cos 9L') (9) 

gives X a lepton number — 1. The interactions in our model arc sufficient to keep the dark matter in thermal 
and chemical equilibrium. 

We define the asymmetries in each field as An^i , An x , and Anx ■ The asymmetries in the fields are (for 
T <M) 

An L , = ^ L T 2 , (10) 



Anx, = -Anx = e~ M ' T (M » T), (11) 

= \^lT 2 (M<T). (12) 

where /i^ is the chemical potential for lepton number. (1l is fixed because the baryon and lepton numbers 
have already equilibriated due to sphaleron processes active prior to the electroweak phase transition. At the 
time of the second stage of the phase transition, the weak gauge boson masses are large enough that sphaleron 
processes are exponentially suppressed. 

These relations are valid as long as the mixing is turned on. At a temperature Tn, bubble nucleation 
occurs and the system moves to the true vacuum via a strongly first-order phase transition. The dynamics of 
the phase transition are discussed in Appendix [B] Therefore, in the moment immediately following the phase 
transition, the particles are still in the same state that they were in prior to the phase transition (according to 
the instantaneous approximation), but their masses are now fixed by the new vacuum. In particular, the mass 
mixing has shut off. The abundances of A, A, and L° arc fixed by the abundances of the mass eigenstates at 
decoupling. In particular, 

An L = cos 2 9Ani J i + sin 2 9An x ,, (13) 
An x = sin 2 6»A?i L / + cos 2 9An x , , (14) 
An x = -An x , (15) 
We now define the total dark matter and baryon numbers (obtained from summing over all charged fields) as 

A x = An x ~An x , (16) 
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A B = An B - An B . 



(17) 



The ratio of X number to An^ is 

A x (1 + cos 2 0)Anx> + sin 2 OAn^i 



An L cos 2 6An L , + sin 2 6 A 



(18) 



A"' 



Once the phase transition completes, there is a conserved global U(1)a"-* symmetry. Since An$ = at the 
time of bubble nucleation, and $ decays to X at late times, the zero temperature value of Ax is identical to 
that immediately after the phase transition occurs. 

Depending on parameters, the ratio of the X and L asymmetries in (|18p can be consistent with both light 
dark matter with a mass of ~ GeV, or with heavier dark matter masses in the range 200 — 500 GeV (other 
models consistent with heavier asymmetric dark matter can be found in [Sl[6j|3]). For the heavy case, the final 
asymmetry ratio is dependent on which terms dominate (1181) . We identify three distinct limits: 



1. Relativistic X. In the relativistic limit, when Tn 3> mx, mxL, there is no thermal suppression in the 
X\ distribution and we have Aul 1 = An X ' ■ Equation (|18[) reduces to 

(19) 



An L 3 

This is consistent with light dark matter, mx ~ 3.3 GeV, which we calculate from (|23p below and by 
imposing Qx = 50b- 

2. Thermally-suppressed X. In the non-relativistic limit mx S> In, then An x , <C An^. In the 
thermal-suppression-dominated limit, An^, / ' An^i 3> tan 2 9 and 



A x „ I2M 3 



This relation is the same as the one appropriate to asymmetry transfer via higher-dimensional operators 
that freeze out when X is non-relativistic [8j. This is the relevant limit when the (j> VEV at Tn or the 
Yukawa coupling yx is small. 

Mixing-angle-suppressed X. The other possible non-relativistic solution is when An X i / An^ <C 
tan 2 9 and the solution is mixing-anglc-dominatcd. In this case, the final ratio of X/L asymmetries is 
independent of An^/ and is determined only by the mixing angle, 

tan 2 9. (21) 



An L 



This is a novel relationship and unique to models with mass mixing. Unlike the thermal-suppression- 
dominatcd limit, this result is only weakly dependent on Tn, since the expression for the <fr VEV (|130D 
is effectively temperature independent at low T. 

We show below that, when the mixing is large (yx ^ 1), both the mixing angle and thermal suppression terms 
in (|18l) are relevant. 

We now relate the baryon number to the X number using the relationship between baryon and lepton 
numbers in the early universe [13] . 

A_b _ 4 6N t + 3 

A L 3 14V f + 9' 1 j 

where N{ is the number of Standard Model flavors. The baryon number to X number in each of the above 



1. 



5 



2. 



= ? (24 ) 



A x 3 tan 2 9 



(25) 



A B 4(1- tan 2 9)' 

2.3 Field content and interactions 

We now describe our model in more detail. It contains: 

• The Standard Model fields, including and most relevant, the Higgs H and left-handed leptons Li. 

• The dark matter fields, which are three vector-like fermions Xi carrying a lepton flavor index. 

• An additional doublet scalar $ with hypercharge +1/2. The mass of $ satisfies to$ > nix so that X is 
stable. 

There exists a Z 2 symmetry, under which $ and X are charged but Standard Model fields arc not. This makes 
X absolutely stable at late times, since ($) = in the true vacuum and the Z 2 symmetry is unbroken. The 
symmetry also excludes the term HG>', which would give $ a tadpole when H gets a VEV. 
The new terms of the Lagrangian are 

C D mxXiXi + y x S-X^ + V(H, *) + h.c, (26) 

where 

V(T = 0) = Ah\H\ 4 - 4 j u 2 |#| 2 + 4fc 2 |$| 4 - 4^ 2 |$| 2 + 4fc 3 |$| 2 |ff| 2 , (27) 

using the notation of [10| that is convenient in the basis of real scalar fields. In addition to ensuring dark 
matter stability, the Z 2 symmetry excludes an term, which would give $ a tadpole when H gets a VEV. 

With the convention in (|27j) . pcrturbativity requires that ki,k 2 < 0.8 and k 3 < 3.1. There are additional 
constraints to avoid hitting a Landau pole below the GUT scale, but this is highly dependent on what other 
fields couple to H and $ (for example, the quartic terms could have large contributions from integrating out 
other weak scale particles) so we do not consider this as a constraint. 

At T = 0, 4> must be stabilized at the origin, while H should have its measured VEV of v = 246 GeV. The 
vacua of the theory are determined by finding the critical points of the potential. To do so, it is simplest to 
move to the basis of real fields. Because U(l) cm remains unbroken, the VEVs can be rotated purely into the 
neutral components of H and $. Furthermore, assuming that the potential is CP-conserving (for simplicity), 
only the real components of $ and H can get VEVs. Therefore, we make the substitution 

(28) 

giving the zero temperature potential 

V(T = 0) = kih 4 - 2njh 2 + k 2 (/) 4 - 2^</) 2 + k 3 (/) 2 h 2 . (29) 

Because we are interested in the evolution of the field values as the universe cools, we also need to include 
thermal contributions to the masses of H and $. These thermal corrections are proportional to T and dominate 
over the tree-level masses for T 3> n. Ultimately, these corrections determine how the phase transition proceeds. 
We can write the finite temperature potential as 

V(T) = k x h 4 - 2fi 2 h 2 + k 2 (/) 4 - 2nl<j> 2 + kz<t> 2 h 2 + -anT 2 h 2 + ^a 2 T 2 <j) 2 , (30) 
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Figure 1: Unshaded regions in the k% — &3 plane give rise to a two stage phase transition and a viable dark 
matter asymmetry through mass mixing for light dark matter. Each parameter point is consistent with dark 
matter of mass 3.3 GcV. The other parameters are held fixed at [i 2 = 54 GeV for both plots, yx — 0.5 (left) 
and yx = 1-7 (right). In the left shaded region (red), the barrier between vacua is too large for rapid bubble 
nucleation and the system remains trapped in the false vacuum. In the right shaded region (also red), the 
4> 7^ direction is not a stable vacuum. In the bottom shaded region (green), the scalar mass is excluded 
by LEP. 



where 

2 e 2 (l + 2cos 2 6%) , 1 2 rqn 
a 1 = 2fc 1 + -fc 3 + + -y t , (31) 

o sin 26% ^ 

,u , 2 ,- , e 2 (l + 2cos 2 g w ) , 1 a f . 

a 2 = 2fc 2 + + r^— + -y x . (32) 

3 sin 26% o 

The light quark and lepton Yukawa couplings are sufficiently small that only the top Yukawa needs to be 
considered for corrections to the h mass. 

We give the details of the vacua and constraints on the parameters in Appendix [A] The first of the two 
most important results is the constraint leading to a two stage phase transition, 



— < — < — 

k 2 Hz a 2 



< % < -> (33) 



which comes from requiring that the <p direction become unstable at the origin before h, while also requiring 
that the Standard Model Higgs vacuum have lower energy than the <p ^ vacuum. The other main requirement 
is that 4> be stabilized to the origin in the true vacuum, 

ml = 2^-^ 2 2 >0. (34) 



2.4 Numerical results 

We present numerical results for a scan of parameters leading to a two stage phase transition. The most relevant 
result for us is the temperature Tn at which bubble nucleation occurs for a given set of parameters. This 
determines which of the three limits (relativistic, thermal-suppression-dominated or mixing-angie-dominated) 
is most relevant and hence the dark matter mass consistent with ^dm = 5f2b ar yon- We present details of the 
calculation of Tn m Appendix [Bl 

The necessary and sufficient conditions for asymmetry transfer through a two stage phase transition are: 

• The constraints (|3"3")) are satisfied. 
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Figure 2: Relationship between the temperature of the first-order phase transition Tn and the <j> quartic 
coupling fc 2 - The other parameters are held fixed at kz = 1.5, yx = 1.7 and \x% = 54 GcV. No first-order phase 
transition takes place in the shaded regions. 



• The bubble nucleation temperature Tn > 0, and the wall velocity is large enough that the transition is 
strongly first-order (see Appendix [Bj. 

• The mass > mx in the final vacuum so that X is stable, where is given in (|34|) . 

There is also a constraint coming from LEP searches for new doublet scalars (specifically slcptons, from which 
there is a bound > 90 GeV [H]). 

The Standard Model Higgs mass parameter, fii, is fixed to be /ii = 42 GeV by a Higgs boson mass of 
120 GeV in the final vacuum, while the quartic coupling ki is fixed by the VEV (h) = 246 GcV. We have a 
four-dimensional parameter space, consisting of the quartic self-coupling k2, the <j>—h mixed quartic coupling 
&3, the (j> mass parameter ^2, and the Yukawa coupling yx- 

The allowed values of /x 2 come from the constraint ([55)1 . The range is ^ 2 = 25 — 100 GeV, which is precisely 
the range we expect if fi\ and /z 2 have a common origin. Since perturbativity limits fc 2 < 0.8, the upper bound 
of fi2 < 100 GeV follows from ([33]) . The lower bound /j, 2 > 25 GeV comes from the fact that the ratio ai/a 2 
is fixed by the Standard Model interactions when yx,k2,k 3 — > 0. The parameter space for a viable two stage 
phase transition is small close to either bound on /i 2 , so we present all of our results with an intermediate, 
representative value of /i 2 = 54 GeV and scan over other parameters. 

The Yukawa coupling is selected to be 0(1). In this case, the symmetric component of X annihilates 
through the Yukawa interaction (see Section [2~5l . 

Wc now present, in turn, the parameters consistent with the three dark matter scenarios discussed above. 

1. Relativistic X. In general, Tn 3> 3.3 GeV, so light 3.3 GeV dark matter is consistent with any choice 
of parameters leading to a two stage phase transition. We show in Figure [1] a region of parameter space 
in the fc 2 — k% plane (with /x 2 = 54 GeV and two values of yx = 0.5, 1.7) with a two stage phase transition 
and light dark matter. 

In the right shaded region, the couplings fc 2 give a large contribution to the thermal mass of <p, stabilizing 
the 4> direction and leading the system to condense directly into the true, h 7^ vacuum. In this case, 
the 4> VEV is never non-zero. For smaller values of fc 2 , the VEV (<fi) sa yLt 2 /-\/fe increases and the barrier 
between vacua becomes larger. In the left shaded region, the barrier is so large that bubble nucleation 
never occurs and the system remains stuck in the false vacuum. The tunneling temperatures are shown 
in Figure [2] 

In the shaded region at the bottom, from (fM]) is smaller than the LEP bound. 

2. Thermally-suppressed X. The relation (|24|) dictates the X masses consistent with a given tunneling 
temperature Tn- We present our results in Figure [3] as a function of the parameters /c 2 , £3, with // 2 = 54 
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Figure 3: Regions in the fc 2 — &3 plane where the final X asymmetry is thermally suppressed. The unshaded 
region is the viable parameter space. Other parameters are held fixed at fj, 2 = 54 GcV and yx = 0.5. Dashed 
contours show representative X masses associated with each point in parameter space (from left to right, 325 
GeV, 375 GeV, 425 GeV). The outer shaded regions to the side (red) do not give rise to a two stage phase 
transition, while the shaded bottom region (green) violates the condition mx < m<j>, rendering X unstable. 



GeV and yx = 0.5. The Yukawa coupling is small enough that the mixing angle contribution to the 
asymmetry is always smaller than the thermal terms. For Yukawa couplings yx > 1, both mixing angle 
and thermal suppressions are relevant and we address this case in point #3. 

The entire unshaded region is consistent with dark matter with a thermally-suppressed asymmetry, and 
the dashed contours illustrate a few representative values of mx between 300 and 450 GcV. The outer 
side shaded regions represent the same constraints as the relativistic case. For small values of k^, 
becomes smaller than the value of mx giving the correct dark matter density, making the dark matter 
unstable. In the inner shaded region on the left, the mixing angle contribution to Ax/Ab dominates 
over the thermally-suppressed terms. 

We find a range of dark matter masses 300 — 500 GeV for the given parameters. When we scan over 
other values of fi 2 , we find a range ~ 250 — 550 GcV, although most of the parameter space lies in the 
300 — 450 range. According to (|33|) . a smaller value of /^ 2 gives a two stage phase transition with smaller 
values of fc 2 , shifting the band in Figure [3] to the left (to smaller values of fc 2 ). Smaller values of yx 
decrease the value of a 2 , allowing larger values of fc 2 and relaxing the upper bound on fc 2 from p3[) . 

This range of masses is the expected range for the two Higgs model since the bubble nucleation temper- 
ature Tn < 2/ii from (|147j) . and this restricts < 550 GcV. For stability, the X mass can never exceed 
the <j) mass, and the largest possible value of m§ with pcrturbative quartic couplings is about 550 GcV. 
The lower bound on masses comes from the fact that the vast majority of points in parameter space 
with a two stage phase transition have tunneling temperature Tn > 35 GeV (when the vacua become 
degenerate at a temperature lower than this, the energy barrier is typically large enough that the system 
remains trapped in the false vacuum). Thus, when the final asymmetry in X is determined by thermal 
suppression, the range of allowed dark matter masses is 250 < mx < 550 GeV. 

3. Mixing-angle-suppressed X. Wc present our results in Figure @] as a function of the parameters A; 2 , 
&3, with /i 2 = 54 GeV and yx = 1.70. Dark matter is consistent with all unshaded regions, with dashed 
contours showing representative values (425 GeV, 450 GeV, 475 GeV) . The outer shaded regions do not 
give rise to a two stage phase transition for the reasons outlined in the relativistic discussion. In the 
lower shaded region (green), the X mass consistent with the asymmetry transfer is much larger than 
m r j >1 making the dark matter unstable. 

1 For yx =0.5 considered in the first two cases, the mixing angle is always too small to give a significant contribution to the 
asymmetry. 
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Figure 4: Regions in the ki — k% plane with yx = 1.7 and {12 = 54 GeV consistent with weak scale dark 
matter. The mixing angle contribution is relevant for all points. The unshaded region is the viable parameter 
space. Dashed contours show representative X masses associated with each parameter (from left to right, 425 
GeV, 450 GeV, 475 GeV). The two inner, unshaded regions show parameters where the mixing angle (left) 
and thermal suppression (right) are respectively dominant. The outer shaded regions (red) do not give rise to 
a two stage phase transition. In the lower shaded region (green), the mx < m-tp condition is violated. 



We mark the regions where the mixing angle and thermal suppression terms in (|18|) dominate. Dark 
matter masses range from 400 — 500 GeV with the parameters chosen. Comparing with Figure El the 
presence of the large mass mixing significantly changes the asymmetry, even in regions where thermal 
terms are largest. This favors heavy dark matter, O(400 GeV), as the large X mass suppresses both the 
mixing angle and thermal contributions and results in the correct dark matter abundance, even in the 
presence of two sources of asymmetry transfer. 

The range of parameters with dominant mixing angle is much more restricted than for relativistic or 
thermally-suppressed dark matter. The reason is that, except for small hi, tunneling occurs shortly after 
the vacua become degenerate and when the <\> VEV is small. This typically leads to a mixing angle that 
is smaller than the Boltzmann suppression factor. 

To summarize, with a two Higgs model with a two stage phase transition, the relationship between the 
baryon and dark matter asymmetries can be determined by thermal effects and/or the mixing angle. With 
large mass mixing {yx ^ 1), the mixing angle contribution to the asymmetry is always important, and the 
ratio of Ax / Ab is modified from the predictions of earlier work where only thermal effects are relevant to the 
dark matter asymmetry [51 [B]. The model is consistent with a wide range of dark matter masses, from 3.3 
GeV in the relativistic case to 300-500 GeV in the thermally-dominated and mixing-anglc-dominatcd regimes. 

Note that for this model the transfer is always in equilibrium and the thermal scattering rate is always 
comparable to the oscillation rate. 

2.5 Annihilation of the symmetric component 

In addition to generating the correct dark matter asymmetry, the coupling must be large enough that the 
symmetric component annihilates efficiently so that the late-time X density is dominated by the asymmetric 
component. Our model provides a built-in mechanism for symmetric X annihilation via the i-channel scalar 
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exchange process process XX^ — > LL^ . The cross section for this process is 



4 2 

_ Vx m x 
167r ml 



(35) 



The requirement that the symmetric number density of X at freeze-out be less than the asymmetric component 
constrains the Yukawa coupling. For masses in the 0(100 GeV) range, this requires x = mx /T{ rcczc ~out ~ 25 
and, 



For values mx ~ m^ ~ 0(500 GeV), this gives the constraint yx > 0.45. This is satisfied for the parameters 
considered above and is true over much of the parameter space consistent with asymmetric dark matter where 
the X asymmetry comes from mass mixing. 

When mx = 3.3 GeV, the suppression in the cross section from (mx/m^) 2 can be very large. When is 
at its lowest bound of 90 GeV, then yx > 0.6. For larger values of m^, then yx > 1 to ensure the annihilation 
of the symmetric component. There may also exist other annihilation channels (for example, if there is a U(l)' 
under which X and <j> are charged) to further deplete the symmetric component. 

2.6 Phenomenology 

While the dark matter in this model is a singlet under the Standard Model gauge group, the model also 
contains a new charged field, the doublet scalar <\>. This leads to a loop-suppressed coupling of X to nucleons, 
as described in [T51 [T] . The direct detection cross section is the same as [5] , 



where m^/yx « 10- 1000 GeV in the two Higgs model (m w 100 - 500 GeV and y x « 0.5 - 10). This 
cross section should be tested by the next generation of detectors. The one-loop interaction with nucleons is 
dominant over direct interactions with electrons in the detector, due to the small mass of the electron and 
therefore the small recoil energy [16] . 

The doublets can be pair produced at the LHC by s-channel weak gauge bosons. The neutral component 
</>° decays to X + u, while the charged components decays to X+£ ± . The decay of ± — > cjp + W ± is suppressed 
by the small mass splitting (induced by electroweak symmetry breaking loops) between components of the 
doublet. 

The charged component is easiest to find. Given the approximate lepton-V flavor symmetry, the signature 
is two oppositely-charged leptons of the same flavor plus missing energy. This is exactly analogous to searches 
for charged left-handed sleptons in the MSSM, with X being analogous to the neutralino LSP. In the case of 
sleptons, assuming no degeneracy with the SLP (m x o <C mj), the LHC will be capable of finding left-handed 
sleptons at the LHC at the 5a level with masses up to 350 GeV with an integrated luminosity of 100 fb _1 at 
14 TeV [T7] . In our case, the doublet masses are in the range 100 — 500 GeV, so there is the possibility of 
detection at the LHC, although discovery will require a long running time. 

This model can be differentiated from supcrsymmetric models because the decay chains should be different. 
Furthermore, this is a more challenging search than gluino or squark searches, so by the time of detection of 
the doublets, we should already know whether or not supersymmetry is present at the weak scale. 

Other possible signatures are the production of two (j>°, which results in a monojet + missing energy 
signature, or the production of </> + </>°, with a signature of one lepton + missing energy. Both are difficult to 
distinguish from the Standard Model background, and we find it unlikely that either is promising. 

Our model is an asymmetric leptophilic dark matter model [6] and can therefore have suppressed indirect 
detection signals. The fact that the late-time dark matter density is dominated by the asymmetric component 
means that there may be insufficient anti-dark-matter in the late universe for there to be an appreciable 
annihilation signal. As shown in earlier works, the symmetric component could, however, be replenished by 




(36) 




(37) 



11 



a small X Majorana mass or there could be some remnant of the symmetric component due to annihilation 
freeze-out [T5] . Therefore, the detection of excesses in cosmic spectra does not rule out this model. As shown in 
[6], the annihilation cross section (|35|) for weak scale dark matter is (c ann v) ss 10~ 24 cm 3 /s, which is consistent 
with the PAMELA and Fermi anomalies for mx = 400 GeV, = 500 GeV and yx = 1 [6]. 

3 Mass Mixing via Planck-Scale-Suppressed Operators 
3.1 Overview 

Our second and more generic scenario exploits large VEVs that we expect in the early universe. Examples of 
fields that are likely to have large, non-zero values include moduli, flat directions in supcrsymmetric theories, 
and Goldstone bosons. At finite temperature and density, the kinetic energies of relativistic fields also have 
sizable background values. Generically, all such fields should couple to any other existing fields through Planck- 
suppressed operators. Since we do not expect such operators to respect any particular global symmetry, these 
operators should lead to the mixing of dark matter and lepton quantum numbers. While theories with moduli 
and flat directions are most likely supcrsymmetric, mixing due to background energy is also present in non- 
supersymmctric theories. 

The key point to take away for all models of mixing by Planck-suppressed operators is that mixing between 
different sectors does occur due to Planck-scale physics (assuming gauge- invariant combinations of fields). 
Therefore, asymmetry transfer between sectors also occurs if the asymmetry generation happens early in the 
universe's history. This is the dominant effect if the dark and Standard Model sectors are very isolated, and 
is still expected to be present as a sub-leading effect if other fields mediate dark-matter-number-violating 
interactions between the two sectors. 

We consider both fermionic and scalar dark matter candidates. Since lower-dimensional operators lead 
to less-suppressed mixing, scalar dark matter proves to be particularly important when the mixing would 
otherwise be too small, such as when induced by thermal kinetic terms and field VEVs, which are temperature- 
suppressed and therefore generically smaller than those for moduli. For operators involving string moduli and 
other flat directions, both fermion and scalar mixing can lead to a large transfer of asymmetry between sectors. 

For simplicity, we consider models in which the dark matter is a Standard Model gauge singlet in order 
to avoid direct detection bounds. This implies that the field mixing with the dark matter is also a Standard 
Model gauge sing leg One option for the field mixing with dark matter is the right-handed neutrino in a 
model with Dirac neutrino masses, which carries lepton number. While wc choose this example to work out 
the concrete details of mass mixing, the general framework is also compatible with Majorana leptogenesis (or 
GUT baryogencsis) in theories with additional lepton-number-carrying gauge singlets. As will be discussed 
later, models with heavy moduli can dilute the dark matter and baryon asymmetries. In the framework of 
Dirac leptogenesis, the dilution can be compensated by resonance enhancement, which is not generic. A more 
general alternative is to have Affleck-Dine leptogenesis (or baryogenesis) along a flat-direction as a composition 
of singlet L(B) and SM leptons (baryons). Alternatively, the moduli can be light enough to be stable over the 
lifetime of the universe. 

We focus on the case where a B — L symmetry is created via Dirac leptogenesis [T5]. In such models, no net 
lepton number is created but an equal and opposite lepton number is sequestered in the LH and RH neutrino 
sectors. The processes that equilibriate the two sectors have rates that are suppressed by neutrino masses and 
so do not come into equilibrium until well after the time of the electroweak phase transition. The asymmetry 
stored in the RH (s)neutrinos can then be transferred to the dark sector via mass mixing terms and provide 
for a connection between the dark matter and baryon asymmetries. 

For our models, we take the field content to be the MSSM as well as the new components outlined below. 
First, there are the dark matter and singlet lepton fields: 

• Three generations of singlet RH neutrino chiral supcrfields Ni . 

• Chiral supcrfields X and X, which comprise the dark matter. 

2 The dark matter can carry a gauge charge as long as the field with which it mixes is also charged. 
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We also assume in this model the existence of fields responsible for Dirac leptogenesis, as outlined in Appendix 



• At least two generations of heavy vector- like doublets and ip, whose decay generates the B — L 
asymmetry. 

• A singlet field x which is necessary for generating an asymmetry in both LH and RH neutrinos. 

• A U(1)jv symmetry that forbids direct coupling between the LH Standard Model lcpton L and the RH 
neutrino N. It can also permit the annihilation of the symmetric component of X if gauged. The gauge 
symmetry is anomaly-free with charges Q(X) = 1, Q(N) = —1, Q(ip) = 1 and Q(x) = 1 with and two 
additional fields, either in the lepton or dark sectors. 

Note that the new matter content is similar to that for models of Majorana leptogenesis (and a dark matter 
sector). But in this case the right-handed neutrinos are light and additional fields are present to generate the 
B — L asymmetry through their decay. 

The tree-level, renormalizablc supcrpotcntial is 

W = M^ip + VuQH u u + y d QH d d + y N NH u ijj + y L L X i>- (38) 

The model possesses a Peccei-Quinn-type symmetry that forbids XX mass terms. These mass scales can be 
generated dynamically at later times, possibly through SUSY breaking. We then expect that the X masses 
(for both the fermions and scalars) are ultimately set by the soft SUSY scale. 

3.2 Mass mixing formalism 

In Section [2] we worked in the mass eigenstate basis. This is appropriate because the mass mixing in the 
two Higgs case is nearly temperature-independent around the decoupling temperatures (particularly when 
compared to the rate of the first-order phase transition) , making the mass cigenbasis temperature independent 
as well. The system also consists of a mixture of relativistic and non-rclativistic states, and the mass basis 
adequately describes the Boltzmann suppression of the heavy states. Finally, the system is fully mixed, meaning 
that the rate of the asymmetry transfer given by {7} is in equilibrium (3> H) and the chemical potentials for 
X and N track their equilibrium values. 

By contrast, the VEVs inducing mass mixing for these models with Planck-suppressed operators arc highly 
tcmpcraturc-dcpcndcnt. In fact, there is no phase transition, with VEVs typically proportional to some 
power of T and turning off gradually. Because the mass mixing changes with time, the mass basis is only 
an instantaneous eigenbasis. To avoid continually changing basis over the course of the time evolution, it is 
simplest to use interaction eigenstates. 

At early times when the Planck-suppressed mass mixing operators are active, the fields are also highly 
relativistic. This is precisely the limit in which the standard neutrino oscillation formalism is valid. This 
section is dedicated to the overview of how distributions of particles with mass mixing evolve over time. A 
good review of this material can be found in [50]. We first present fermion mixing and then discuss scalars. 

We first review the story with pure coherent oscillation, which is directly analogous to neutrino oscillation. 
For Weyl fermions ip a , a = 1, . . . ,N, the mass part of the Hamiltonian is 

"Hinass = TjMaP^alpp + h-C (39) 

We assume no mass degeneracies in order to simplify the results. We use Greek letters to label flavor eigenstates 
and Latin letters to label mass eigenstates. The mass eigenvalues Mf are found by diagonalizing the matrix 
M'M. The rotation matrix U that diagonalizes M satisfies 

Y,Ul(M^M) a ^=Mf6 ir (40) 
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The relationship between the mass basis and the flavor basis is 



\fl> a )=Y^U ai \i> i ). (41) 

i 

During a particle's propagation between interaction points, it can be described by a wavefunction in 
standard quantum mechanics. After a time t, the probability of oscillation from flavor state a to ft is given 
by V a ^rp(t) = |(^(i)|^ Q (0))| 2 . Expanding the states in the mass basis and performing space-time evolution 
using the operator e - li "+*P- x gives, in the relativistic limit, 



x , (Ml- Mf ^ 



(42) 



where a ^ ft. The case where a = ft is simply found by subtracting the probability of oscillating into anything 
else. 

For the special case of two fields, we can write U in terms of a single mixing angle, 

U=( C0S6 S[U0 ). (43) 
V — sin 8 cos 8 I 

The oscillation probability then reduces to the familiar form 

P a -+p{t) = sin 2 28 sin 2 f 2 — 1 1 j . (44) 

We compute the exact rate at which the asymmetry is transferred between sectors by mass mixing using the 
density matrix formalism for fermions and scalars. Before doing so, we present a good, approximate formula 
for the rate of asymmetry transfer that agrees with the density matrix results in two important limits discussed 
in Sections 13.2.11 and 13.2.21 This formula will be useful in understanding how the asymmetry rate depends on 
the mass eigenvalues and mixing angle. We derived all of our numerical results in Section [3] using the density 
matrix formalism. 

In an environment where frequent scatterings occur, as is typical in the early universe, oscillation is inter- 
rupted by interactions with other states. This decoheres the system and can be thought of as collapsing the ip 
wavefunction into one of its interaction eigenstates. Interactions are occurring at a total average rate To, 

r = ^(r Q + i», (45) 

where r a is the total rate with which species a interacts with the thermal background. The mean free path is 
1/IV The rate that species a changes into species ft is going to be given by this interaction rate, multiplied 
by the probability of oscillation: 

T a ^ = T sin 2 28 sin 2 ^ = T sin 2 28 sin 2 (jL^j , (46) 

where e = (Mf — M\)j2E. Equation (|46p is easy to understand and is a good approximation to the actual 
evolution of various species in a thermal background and we use it frequently. If T a ^p > H, the mixing is in 
equilibrium. Otherwise, we can obtain an order-of-magnitude estimate of the amount of transfer from species 
a to ft by the quantity Ta^p/H evaluated at the time the asymmetry is generated in one sector. 

The a — > ft transfer can proceed in two different ways associated with different limits of (|46)) . Which regime 
the system is in is independent of whether or not the mixing is in equilibrium. When e/To S> 1, the oscillation 
length is much shorter than the mean free path and so a larger amount of a particles are transferred per 
collision. As a result, fewer collisions are necessary to bring the system into equilibrium. The opposite limit, 
e/To 1, occurs when collisions occur very rapidly compared to the oscillation time scale and few a particles 
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are converted to /? particles per collision. Nevertheless, the transfer process can still be in equilibrium if the 
mixing angle is large. We encounter examples of both limiting cases in Sections 13.31 and 13.41 

The evolution of the different flavor states, taking into account oscillations and coherence-destroying ther- 
mal interactions, can be performed in a unified framework using the density matrix formalism \12\ \21\ 122] . 
The system begins in some state po, where the diagonal components give the iniital relative abundances of 
each particle type and the off-diagonal components represent any coherences present in the initial state. The 
time evolution of the density matrix is given by the equation 

ip=[H^,p]-i{H^,p}, (47) 

where is the real part of the Hamiltonian calculated to leading order and includes masses and interactions 
with the thermal background, while consists of the imaginary part of the background potential consisting 
of absorptive interactions with the background that cause the system to decohere. Integrating (|47p yields the 
final abundance of each species. In Sections 13.2.11 and 13.2.21 wc give (|4T|) in component form for scalar and 
fcrmion mixing and show that, in the limit of time-independent temperature and mixing, the result reduces 
to 0BJ. 

Wc arc interested in determining the evolution of an asymmetry that is stored entirely in one sector and 
subsequently migrates to a different sector through oscillations. Since the evolution equations arc linear, it is 
appropriate to consider only the asymmetric components. The initial density matrix at the time the asymmetry 
is generated and that we use in our analysis is po = diag(l, 0, 0, ... , 0), where the 1 is the diagonal component 
of the species housing the asymmetry. 

At finite temperature and density, particularly temperatures where T > M iy finite temperature contribu- 
tions to the background potential must be included in the mass mixing formulae (we assume that both sectors 
are in kinetic equilibrium so that there is a common temperature between them) . Considering again the 2x2 
case, the effective masses from thermal corrections for both fermions and scalars in the interaction basis are 

|M(T)| 2 - |Af(0)| 2 = ( A? q T2 a 2° t2 ^ , (48) 

where Ai and A2 are effective couplings parameterizing all interactions with the background plasm^f). In 
Section [3.3] we relate Ai and A2 to the couplings of specific models. 

Using E w 3T at finite temperature and density in the ultra-relativistic limit, we get the Hamiltonian 

H afi =3TS aP + i \J} aP . (49) 
bl 

In the case that there exist common interactions between the sectors that arc in thermal equilibrium, T is the 
same for all fields. The mixing angle and energy eigenvalue splitting arc unchanged when a term proportional 
to the identity is subtracted, so we can subtract the T 6 a p term from (j49]) . 
We now derive results for the particular cases of scalars and fermions. 

3.2.1 Scalars 

Let the tree- level, diagonal elements of My be denoted by pi. Then, the total mass-squared matrix is 



, u 2 + A 2 T 2 M 2 

' M < T >'=I ui + < 50 ' 



The mass eigenvalues and mixing angle are 



2m 2 ± = \\T + A 2 T + p\ + p\ ± J (A^T 2 + p\ — p\) 2 + 4M* 2 , (51) 



3 For fermions, finite temperature corrections do not give rise to a mass term in the Lagrangian due to chiral symmetry, but 
they do give rise to an effective potential that can be considered as a contribution to the M matrix 
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sin( 



v/2Mf 2 



(52) 



where A|i = A 2 . — A 2 . 

The Hamiltonian JH) is 



and 




p= PU P12 

The density matrix evolution equations (|47j) in component form are |22j 



(53) 



(54) 
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-iM? 2 



In this case, flS]) gives 
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-iM 2 2 



A^T 2 



1 



r = -(Ti + r a ) 



iM* 2 



-2T T - i [nl 



Pi 



AiiT 2 ] J 





(Pll\ 




P22 




P12 


) 


\P21J 



(55) 
(56) 



Note that the previous formula applies also in the limit that M\ 2 and T are independent of time, in which 
case we have a linear system of first-order differential equations. The solutions are exponential functions, 
with the coefficients in the exponential being given by the eigenvalues of the matrix in (|55|) . There is one 
zero eigenvalue, representing the equilibrium distribution, and the approach to chemical equilibrium (i.e. the 
asymmetries being equal in both fields) is given by the smallest of the eigenvalues, denoted by 7 s iow Defining 



6T 

the solutions can be found analytically in the limits To ^ e and To S> e. The solutions are 



(57) 



7 s i ow = 2r shr 2 9 cos z 6, T Q < e, 

7siow = sin 2 6 cos 2 6—, T a > e, 
1 



(58) 
(59) 



which reproduces (|46|) in the same limits. This confirms our physical picture of transfer via mass mixing and 
interactions with the background, and also checks equation (|46[) . 

Since all of the parameters in our theory are actually time-dependent, the system never reaches chemical 
equilibrium if 7 s i ow < H . There is a significant parameter space for the models in Sections l3.3l and l3T4"l for which 
this is true. In this case, the amount of asymmetry transferred between sectors is approximately j s i ow /H. 



3.2.2 Fermions 

Suppose we have three Weyl fermions with a Dirac mass /x between ipi and ij)2 and a Dirac mass mixing M13 
between ipi and ^3. The mass-squared matrix is 



\M{Tf 



M 2 3 





\\T 2 







A 2 T 2 



fiM- 



13 





fiMia 
M 2 3 + A|T 2 



(60) 
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Unlike in the scalar case, the tree-level mass p appears in the off-diagonal terms. The matrix is block diagonal 
and only the lower 2x2 sector is mixed. The mass eigenvalues and mixing angle are 



2m| 



sin 9 



fj, 



Mf 3 + (A 2 + X 2 3 )T 2 ± V^ 4 + (M? 3 + \l 2 T*y + 2 M 2 (M 2 3 - A§ 2 T 2 ), 



y/2/J.M. 



4 13 



(61) 
(62) 



where 

The Hamiltonian is 



A = - A 2 2 T 2 ) 2 + 2M^ 3 (^ + A 2 2 T 2 ). 



H = ±- 

6T 



( p 2 + M 2 3 + \\T 2 







\ 



V 



(63) 
(64) 



(65) 





p? + X Z 2 T Z pM 13 

fj,M 13 Mf 3 + a|t 2 y 

The density matrix is 

^Pn P12 P13 
P = P21 P22 P23 

\P31 P32 P33 

Because H is block diagonal, the density matrix evolution equations involving ipi decouple from those inde 
pendent of ipi, giving three independent systems of equations, 

/ ipM 13 

-ipM 13 

ipM 13 -ipM 13 -6r T + i(n 2 - M 2 3 - „ 32 
\ -i[iM 13 ipM 13 



( P22 ^ 




P33 


1 


P23 


" 6T 


\f>32 j 




( P12 


| = — 


\P13J 


1 6T 


fill = 0. 



A 2 T 2 ) 



-ipMi 3 
ipM l3 


-eroT-i^-M^-Al^y 
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( P22\ 




P33 




P23 


/ 


\P32J 



1 / -6T' T - iM{ 3 ipM 
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1:; 



-6r^'T - i^i 2 



P12 

Pl3 



(66) 



The equations for p 2 \ and p 3 \ are found by taking the conjugate of those for p\ 2 and p\ 3 . Because there are 
now three particle types, we have three rates and 



r 



1 



(r 2 + r, 



To = ^(ri + r,), 



1 n 



2( r i+ r 3) 



(67) 
(68) 
(69) 



Since pn is conserved, we can only have transfer between particles of species 2 and 3, and therefore only the 
first system of equations in is relevant for our purposes. 

As in the scalar case, the analytic solution to (|6"6")l in the case of time-independent T and M 13 is 



7 s i ow = 2r sin 2 9 cos 2 9, T < e, 

2 2 e 2 
Tsiow = sin 6* cos 0— , T > e, 

i 



(70) 
(71) 



which reproduces (l4"6"l) . 
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3.3 Moduli-induced mass mixing 
3.3.1 Overview 

One class of fields that are expected to have large background values at early times are moduli fields; examples 
include string moduli or Polonyi fields for supersymmctry (SUSY) breaking. As their potentials are exactly flat 
at the perturbative level in SUSY, moduli can be found at finite temperature far from their zero-temperature 
minima. Such fields can have VEVs as large as Mp\ following inflation. The moduli subsequently roll down to 
their true minima, which would be established by SUSY-brcaking terms in the potential. During this period 
of rolling, the VEVs can induce large mass mixing terms between other fields in the theory. 

Note that moduli decay can create entropy that would dilute both lcpton and dark matter numbers. We 
assume that either decay occurs very late (as with light moduli) with Dirac leptogenesis or that a more efficient 
means of baryogenesis (such as Affleck-Dine baryogenesis) applies. Dirac leptogenesis can be efficient as well 
with resonance enhancement, but it is less generic. We show that both light and heavy moduli are possible in 
Section [3.3.51 For clarity, we use bold- face when discussing the cases of heavy and light moduli. 

We consider the dark matter to be a vector- like pair of superfields X, X with tree-level mass /itree- Mass 
mixing arises from Planck-suppressed higher-dimensional operators. Mass mixing between X and N generically 
occurs for both fermions and their scalar counterparts. 

For heavy moduli, mixing between fermion components and a modulus <f>, mixing can come from super- 
potential terms such as 

cd 2 

AW — -j—XN + h.c, (72) 
Mpi 

which gives a mixing 

m X N = -T7-- (73) 
Mpi 

Since the Hamiltonian contains the squared mass matrices, the relevant mixing term is not mxN but the 
product 

2 CfJ,tree4 > 

"Vix = Mtree m X N = JTT ■ (74) 

Mpi 

The dark matter is stable if /Xtroo < 2m^. 

Superpotential terms also generate mixing between the scalar components of X and N. In the presence of 
a spurion field S with SUSY-breaking .F-term, an allowed operator is 

c S6 2 

AW = ^—j-XN + h.c, (75) 
If Ts is responsible for SUSY breaking in the MSSM, the mixing is 

2 C s m 3/2 4> 2 ,- a , 

m XN = 77 • (76) 

Mpi 

We compare the mixing term in the squared mass matrix (|76[) with the fermion case (|74[) . The ratio of fermion 
to scalar mixing is /^troo/"T-3/2- 

The late-time dark matter particle can be either scalar or fermion, depending on the details of the model. It 
is, however, irrelevant at late times whether the dominant mixing was through scalars or fermions, since inter- 
actions within the X sector (most importantly, those responsible for annihilation of the symmetric component) 
cquilibriate the asymmetry among the scalar and fermion components. 

We expect that higher-dimensional couplings to moduli fields also contribute to XX mass terms. For the 
fermion case, 

AW= ^-II + h.c, (77) 
Mpi 
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This operator doesn't violate dark matter or lepton numbers, but is relevant to the mixing angle. The resulting 
diagonal mass is 

, CX ^ - i lno\ 

= Mtrcc + -77 = Mtrco + « m X N- (78) 

Mpi 

Since both XX and XN get masses from the moduli VEV, we express everything in terms of mxN and the 
proportionality constant n = m xx /mxN, which we expect to be 0(1) since both masses arise from the same 
moduli fields and operators of the same dimension. For scalars, the corresponding diagonal mass term is 

2 _ 2 , KC s m 3/2 (t) 2 

scalar — Mtrcc ' 77 ■ \ IJ ) 

Mpi 



In the scalar case, the diagonal mass does not contribute to the off-diagonal component of the mass-squared 

' scalar ' 



matrix, rn XN . We show in Section 13.3.31 that the asymmetry transfer rate does not depend on fi x s 



For light moduli, the condition mx < 2ro^, is not satisfied. The fermion decay rate is. 

r*^, (80) 

and the X lifetime would be much below the bound of r > 10 26 seconds [23] . For this model, we assume a Z3 
symmetry that only allows higher-dimensional operators such as 



C0 3 

^ w light = 

which gives a mixing 



AWught = 7^M + h.c, (81) 



M 2 



.3 



C(f> 

The diagonal component of the mass-squared matrix is 



m X N = jfi-- (82) 



2 cfi X (f> ,„„, 
"W = 71 t2 ■ ( 83 ) 



The corresponding superpotential term for scalar mixing is 



AW = -jrg-XN + h.c, (84) 



which gives a mixing term in the mass-squared matrix 



2 _ cm 3 / 2 (l> 

lXN ~^ 2 



.3 



(85) 



Pi 

As with heavy moduli, the ratio of fermion to scalar mixing is jJLx/fns/2- The decay width is 

I* « (86) 

safely satisfying the constraints for m,x ~ GcV-TcV. The mixing term for scalars is the SUSY-brcaking spurion 
S coupled to (|ST]|. 

We present a toy model for the potential of moduli of any mass fM\ : 

V = (m 2 a 2 H 2 )\(f>\ 2 + ^2 {ml + b 2 H 2 )\cj>\ 4 , (87) 
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and H 2 mass terms arise from SUSY breaking due to MSSM SUSY 
breaking terms and finite-temperature effects, respectively. At early times, the minimum of this potential lies 
at 

ja 2 H(t) 2 -m 2 

W=^W , 2 J^ , J - (88) 



b 2 H(ty 



m 



At the critical time t c , when H(t c ) = m^/a, the minimum of the potential disappears. Because of damping 
in the equations of motion, <f> is unable to efficiently track the minimum of the potential near t c and is left 
at some finite value. After the critical time, it begins to oscillate with a power-law-damped envelope, rolling 
toward the true minimum of <f> = 0. Its equation of motion is 

4> + 3H<P + 2(77^ _ a 2 H 2 )<f> + -jL-(m% + b 2 H 2 )<f> 3 = (89) 

M pl 

To obtain the asymmetry transfer rate, we must also specify the interactions in the theory. The simplest 
model that allows for the annihilation of the symmetric component of A is a U(l)' gauge interaction with 
coupling g under which A and N are oppositely charged. In the most minimal scenario, N also has Yukawa 
interactions (with coupling y) from leptogcnesis. We then have, in the notation of section [ 



A| = 5 2 /8, (90) 
X 2 N = g 2 /8 + y 2 /8, (91) 
>?nx = 2/78- (92) 

The Yukawa coupling y depends on the model of leptogenesis. For Dirac leptogenesis, the effective coupling 
is y = M&x(yl,yjj)T 2 cp /M^ with T lcp < M^. For Affleck-Dine leptogenesis, y = (T\ ep /M P \) n , where n is the 
dimension of the operator lifting the flat direction. 

Because of its importance, we repeat the analytic estimate of the transfer rate P6l) : 



r n->x = Tq sin 29 sin 



2n„ • 2 



e 

\2T~o 



Three parameters influence the transfer rate: the mass eigenstate energy splitting e = (rn\ — m_)/6T, the 
mixing angle 9 (both of which are explicitly given in eqs. (|6H62j) ). and the thermal scattering rate Tq w 
8^r(ff 4 + y 2 )T ■ The asymmetry transfer rate is reduced for small mixing angles and small mass splittings. 

Note that in the limit where e ^> To, the oscillation is rapid and sin 2 f^T^) — > 1/2 takes on its average value. 
There are constraints associated with the couplings and scales in theory. They are: 

• Leptogenesis should occur after reheating to avoid diluting the lepton asymmetry with the massive 
quantity of entropy released by inflaton decay. 

• In order to ensure efficient depletion of the symmetric component of A, the A annihilation coupling 
should be g - 0(1). 

• y <C 1 for the heavy field in leptogenesis to decay out of equilibrium. 

• T > mxN to avoid thermal suppression of A and N fields. 



3.3.2 Evolution of </> 

The evolution of 4> ([89)) depends on its mass. We present the evolution of <j> for both heavy moduli (m^ > 50 
TeV) and light moduli (ra^ < 10 MeV). We justify these mass ranges in Section [3.3.51 

We first consider heavy moduli, <fr oscillation begins when H ps m^. For > 50 TeV and Trh ^ 10 12 
GeV, oscillation will begin prior to inflaton decay. In between the time t c when oscillation begins and reheating, 
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the cosmic background density is matter dominated by inflaton and moduli fields that scale as massive matter. 
This gives an approximate solution to the 4> evolution: 



sin (m^t) 

4>~4>o—, t^-, (93) 

where 0o ~ Mpi. 

At the critical time t c , H w and (4>) = Mpi, so the <f> and inflaton fields have comparable energy 
densities. Before the inflaton decays, both p^ and p- ln f scale identically with T, so their densities are also 
comparable at the time of reheating. Setting p^ = pj n f at the time of reheating establishes the relationship 
between temperature and time in the subsequent 0-dominatcd stage of the universe's evolution. Assuming the 
inflaton decays entirely to radiation, equating the inflaton and moduli densities gives 

^S*T£h « m^(i RH ) 2 . (94) 

The background remains matter dominated and the relation (|93[) stays valid until the time of moduli decay, 
and during this epoch the temperature-time relation is 

T = T RH (Wi) 2/3 . (95) 
The solution to the evolution of </> for heavy moduli is 



The mass mixing follows from ([73)) and (|76|) . Asymmetry transfer from N to X begins only at the time 
at which leptogenesis occurs, t\ cp . Because mxN turns off as T 3 , the dark matter asymmetry is essentially 
established at ti op . The mass mixing values mxN clearly depend on the model input parameters a, b, c, and 

We now turn to light moduli. In this case, oscillation must begin before H — to avoid having <f> 
overdose the universe. One possibility is that the <f> negative mass-squared predominantly comes from inflaton 
couplings that decay at reheating (see Scction [3.3.5p for details. Oscillation then begins at the time of reheating. 

For all values of T h h and light satisfying the overclosure bound (|115p . oscillation begins only after 
reheating. As a result, the universe is radiation-dominated during the time of leptogenesis and mass-mixing, 
in contrast to the matter-dominated scenario for heavy moduli. In a radiation-dominated universe, T and t 
are related by [55] 

. - ^ (97) 



and the 6 evolution is 



The mass mixing is obtained from (|82[) for fcrmions and (|85[) for scalars. 
3.3.3 Analytic approximations 

We derive approximations that directly link the parameters from moduli mixing to the mixing angle and mass 
eigenvalue splittings in Section 13.21 There are three relevant limits for the mass eigenvalue splitting e and 
sin6>. These are mxN S> fJ-x, ttixn "C px, and mxN ~ f^x- 
The results for fcrmions follow from equations (|6lT) and ([62]) : 



e=^|, sin0 = ™, p x -»mx N ; (99) 
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2 

XN 



, S in6 = fx x <^m XN ; (100) 



^f, sin0=^, ^x=m XJV . (101) 



The overall thermal interaction rate is To « <? 4 T/87r 3 . When the gauge coupling is 0(1) (as needed for 
annihilation of the symmetric component of dark matter), ttixn, fJ-x < T in order to avoid thermal suppression. 
This means that typically e < To and the asymmetry transfer rate in each limit (|46[) reduces to the same 
quantity, 

_ 4 sin 2 cos 2 6>e 2 _ /x| _ ^</> 4 

2r ~ 18T 2 r ~ 3T 2 T M^ [ ! 

For k = 1 (i.e. including higher-dimensional corrections to the diagonal mass) and ttlxn 3* Mtroej this rate 
> H and the system is in chemical equilibrium for mxN ~ T, while it is out of chemical equilibrium for lower 
mixings. In this limit, 

M 33* iw*," (103) 

Since </> is proportional to (T/Trh) 3 , if reheating and leptogenesis occur around the same scale, the asymmetry 
transfer is in equilibrium and chemical equilibrium is reached. If, however, there is a large hierarchy between 
Trh and T\ ep , the transfer rate is suppressed. 

In the case of scalars, from equations ([51"]) and (|52|) : 

2 

e=^, sinfl=^, Mx»m X jv; ( 104 ) 

o 

Tfl 

e =~ff-> s[ne= 72> Mx<mxjv; (105) 
e=^K sine=^=, Hx = m XN . (106) 

Oi V5-V5 

For similar reasons to the fermion case, when ttixn, (J-x < T, the asymmetry transfer rate (|46[) in all three 
limits is approximately 

4sin 2 flcos 2 g £ 2 ^ m 4 ^ ^ m§ /2 4 
^ X 2r ~ 18T 2 r ~ 18T 2 T M^ { ' 

Once again, when leptogenesis occurs close to the scale of reheating, ttixn is large and the system is in chemical 
equilibrium. At lower leptogenesis scales and smaller values of mxN, the system is out of chemical equilibrium. 
We can compare the transfer rate from mixing between fcrmions and scalars: 

2 

fiV-^X, scalar ^3/2 (108) 



i 2 
N^X, fermion Mx 



Therefore, scalars give the dominant contribution to mixing if m 3 / 2 3> fJ-x, and fermions give the dominant 
contribution if fix 3> rn 3 / 2 . Of course, the heavier of the two state will decay to the lighter state, which is the 
dark matter. 



3.3.4 Numerical results 

To confirm our approximate results, we numerically solve the density matrix evolution equations to determine 
the magnitude of the TV — >• X asymmetry transfer. We summarize the results below. The dark matter density 
is firjM ~ (5 — 6)r^B- We choose parameters so that IJdm = 5f2s- 

We begin with heavy moduli. <j> oscillation begins before reheating for > 50 TeV, and 4> is damped 
from its initial value of Mj>\ by the time of reheating. For Trh ~ 10 10 GeV, reheating occurs only shortly 



22 



after oscillation begins and 4> is n °t much damped from Mpj. As a result, fermion mass mixing is > T and 
fcrmionic modes are suppressed, allowing the scalar mixing to dominate the asymmetry transfer. 

With smaller reheating temperatures, T ~ 10 8 TeV, <j> is sufficiently damped at reheating that fermionic 
modes are not thermally suppressed. In this case, either scalar or fermion mixing can dominate, depending on 
the ratio of nx/m 3 f 2 , according to ()108[) . 

We present our results in terms of benchmark points for a high reheat temperature of Trh = 10 10 GeV 
and a lower reheat temperature of Trh = 10 s GeV. We use — 50 TeV for both cases. 

For Trh = 10 8 GeV and all leptogenesis scales T ep < Tan, fermion mixing nixN < T and so we must 
consider both scalar and fermion mixing. With a = b = 1, c = k = 2, some benchmark points are 

1- T cp = 10 8 GeV and mx = 5 GeV. The dominant mixing is from fermions because 
Hx ~ 5 x 10 5 GeV > m 3/2 5 x 10 4 GeV. 

2. T C p = 10 7 GeV and mx = 90 GeV. The dominant mixing is from scalars because 
Hx ~ 90 GeV < m 3/2 . 

The transfer rate is fast enough to bring X and N into chemical equilibrium for point #1 above, consistent 
with ~ 5 GeV dark matter. The transferred asymmetry is suppressed for point #2, consistent with heavier 
dark matter. Heavier masses > 100 GeV are consistent with later leptogenesis times. 
For Trh = 10 10 GeV, with a = 6 = 1 and c = k = 0.1: 

• For leptogenesis scales 5 x 10 8 GeV < Ti cp < Trh, fermion mass mixing is > T and fermionic modes are 
suppressed. Mixing from scalars gives: 

T lcp = 10 8 - 10 10 GeV and m x = 5 GeV. 

• Below Top = 5 x 10 8 GeV: 

1. Ti e p = 5 x 10 8 GeV and mx ~ 5 GeV. Mixing is dominated by fermions because 
Hx ~ 3 x 10 8 GeV > m 3/2 . 

2. T op = 5 x 10 7 GeV and mx ~ 5 GeV. Mixing is dominated by fermions because 
Hx ~3x 10 5 GeV > m 3/2 . 

Scalar mixing dominates for Ti ep < 10 GeV and asymmetry transfer is suppressed, consistent with X 
masses up to 100 TeV depending on the leptogenesis time. Note that, even though the gap between 
the reheat temperature and leptogenesis scale is wider in this example than for the low reheat case, the 
4> VEV (and hence mass mixing) is still large enough to keep the system in chemical equilibrium at 
Tcp = 5x 10 7 GeV. The reason is that, with a higher temperature, the reheat time is also earlier and so 
the time of leptogenesis is earlier than it would be with a lower reheat temperature. Since <fi ~ 1/t, this 
leads to a larger mixing. 

To summarize our findings for heavy moduli, asymmetry transfer from higher-dimensional moduli couplings 
to XN is consistent with dark matter masses ranging from 5 GeV to the weak scale over a range of reheat 
temperatures and leptogenesis scales. We generally expect scalar mixing to dominate when mxN > T and 
the fermion modes are suppressed. With lower T ep , fermion mixing dominates over scalar mixing when mxN 
is just below T, while scalar mixing dominates once again at low T ep because mxN is small and the scalar 
mixing has an ra 3 / 2 enhancement. 

Moving on to light moduli, <j> = Mpi at the time of reheating. Therefore, if leptogenesis occurs at or 
near the reheat temperature, the fermion X and N are thermally suppressed and the dominant asymmetry 
transfer happens between scalars. The scalar mixing is small because m 3 / 2 ~ m^ is small in this case, and the 
resulting transfer is out of equilibrium. In fact, for the points shown below, mixing from scalars alone cannot 
transfer enough asymmetry to be consistent with a dark matter mass below 100 TeV, which is the upper limit 
for fields whose symmetric components can be annihilated by perturbative couplings. 

We present benchmark points for reheat temperatures Trh = 10 9 GeV and Trh = 10 11 GeV. We begin 
with Trh = 10 9 GeV and m<p = 3 keV, which is the largest moduli mass consistent with decay and overclosure 
constraints. 
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• For leptogenesis scales 2 x 10 6 GeV < Ti cp < Trh, fermion mass mixing is > T and fermionic modes are 
suppressed. Mixing from scalars is highly out-of-equilibrium and transfers no appreciable asymmetry. 

• For leptogenesis scales T cp < 2 x 10 6 GeV, fermion mixing is dominant. The asymmetry transferred to 
X drops with 2I ep . Some benchmark points are 

1. T lop = 2 x 10 6 GeV and m x = 5 GeV. 

2. Ticp = 5 x 10 5 GeV and m x = 5 GeV. 

3. T lcp = 10 5 GeV and m x = 160 GeV. 

For the first two benchmark points, mixing is large and chemical equilibrium is attained as outlined in the 
previous paragraph. The transfer rate is suppressed for lower T\ ep and so Ax "C Ab for the last case. Higher 
X masses are consistent with lower T ep , but in practice, it is very difficult to achieve leptogenesis at such a 
low scale. 

For a high reheat temperature of Trh = 10 11 GeV, with a = b = c = n= \. As before, we use = 3 
kcV. 

• For leptogenesis scales 8 x 10 8 GeV < Ti cp < Trh, fermion mass mixing is > T and fermionic modes are 
suppressed. Mixing from scalars is highly out-of-cquilibrium and transfers no appreciable asymmetry. 

• For leptogenesis scales T cp < 8 x 10 8 GeV, fermion mixing is dominant. The asymmetry drops with 
Ti C p. Some benchmark points are 

1. Top = 7 x 10 8 GeV and m x ~ 5 GeV. 

2. Tcp = 5 x 10 8 GeV and m x ~ 5 GeV. 

3. Tcp = 10 8 GeV and m x = 2.2 TcV. 

4. Tcp = 3 x 10 7 GeV and m x = 80 TcV. 

For the first two points, T^^ x > H and the system is in chemical equilibrium, leading to n x ub and 
m x ~ 5 GeV. For the latter two, the transfer rate is suppressed and n x <C ?ib, consistent with X above the 
weak scale. 

There are two principal differences between the results for light and heavy moduli. Since the modulus 
mass is tied to m 3 / 2 , the ratio of soft scales between light and heavy moduli is ~ keV/TeV ~ 10~ 9 . As a 
result, scalar mixing is much more highly suppressed with light moduli. The second difference is that, for light 
moduli, <fi only begins oscillation at reheating, whereas for heavy moduli oscillation begins before reheating. 
This means that the <p VEV is larger for light moduli, leading to larger fermion mixing. 

For both light and heavy moduli in the limit Tq <C H, the particles do not undergo any scatterings with the 
thermal background in a Hubble volume. This means that X and N remain in a coherent state given by the 
oscillation between the two states. The asymmetry transferred in this case is dependent only on the mixing 
angle. The requirement that Tq <^ H is, however, at odds with the condition that the symmetric component 
of the dark matter be efficiently depleted. Therefore, the scenario of oscillation with no thermal scatterings is 
not a viable asymmetric dark matter model. 

In summary, we have presented a number of scenarios where mixing between dark matter and lepton- 
numbcr-carrying fields is induced by moduli VEVs. Depending on the reheat temperature and scale of lepto- 
genesis, the mixing can be very rapid, bringing the two sectors into chemical equilibrium and transferring an 
X asymmetry consistent with GcV-scale dark matter. Alternatively, the mixing can be suppressed, leading to 
dark matter at the weak scale or heavier. The mixing suppression can be very large, but the largest possible 
dark matter masses are ~ 100 TcV, above which point the symmetric component of X cannot be annihilated 
by perturbative couplings. 
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3.3.5 Moduli cosmology 

Moduli are gravitationally coupled and therefore long-lived, with characteristic decay widths 

(109) 



As a result, moduli tend to dominate the energy density of the universe until the time of their decay. The 
entropy accompanying moduli decay can greatly dilute any existing relics, including the baryon asymmetry. 
Furthermore, if moduli decay occurs after Big-Bang Nucleosynthesis (BBN), the resulting entropy generation 
alters the predictions for the light element abundances and such a scenario is in conflict with observations. 
There are two possibilities: either moduli decay reheats the universe to above the scale of BBN (about 5 McV), 
or they are stable to the present day. These possibilities require heavy moduli and light moduli, respectively. 

Assuming that the energy of the universe is moduli-dominated for the case of heavy moduli, the reheat 
temperature T\ from moduli decay is given by 

/3(WU(t = i/r ) 2 \ 1/4 

T\ = ^ . (110) 

\ * 9* J 

Since T\ > Tbbn ~ 5 MeV, we determine that > 50 TeV. 
The temperature immediately before <f> decay is 

To =T RH (iR H r>) 2/3 . (Ill) 

The generated entropy (and hence the dilution of particle relics) is 

so J 

For 7?i0 = 100 TeV and Trh ~ 10 9 GeV, the dilution factor is O(10 10 ). This requires a primordial baryon 
asymmetry of hb/s ~ 1 to give the correct abundances at late times. 

One method of efficient baryogenesis is Affleck-Dine leptogenesis [26] [24] along a flat direction with non-zero 
N and Standard Model lepton number. If the entropy generated from the decay of the flat direction is larger 
than that from the infiaton, then Ub/s can be of order unity |24j . The baryon asymmetry is 

rri / \ 2/(n-2) 

n B ^ ?ib Trh / m FD \ 
s uft) mpY) V Mpj / 

where ns/riFD ~ i n the standard Affleck-Dine set-up. For mpD ~ TeV, Trh ~ 10 9 GeV, and n > 7, 

then ub/s ~ 0(1) and the dilution from moduli decay leads to the correct baryon asymmetry at late times. 
Examples of a flat direction lifted at dimension-8 and invariant under the Standard Model and U(l)' gauge 
groups include: </>p D = NLLE and NUDD, where we impose i?-parity to forbid a lifting term of dimcnsion-4. 
A detailed account of Afflcck-Dinc leptogenesis (baryogenesis) and moduli decay can be found in [57]. Other, 
highly efficient, non-standard leptogenesis mechanisms (such as resonant leptogenesis [55]) are also possible, 
but less generic. 

If the moduli are instead light, they can be stable on cosmological scales. From cosmic ray data, the 
constraint on the dark matter lifetime is t > 10 26 seconds [53]. To be conservative, we require that the moduli 
also have a longer lifetime than this. Using the moduli decay rate (|109p . long-lived moduli have masses < 7 
MeV. 

There are additional constraints on moduli masses from the requirement that moduli do not overdose the 
universe. If the initial </> VEV is Mp\ and oscillation begins at a temperature T osc , then <j> energy density to 
entropy ratio is 

T * Tf^- (114) 
* y * 1 osc 
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To avoid overdosing the universe, and so that moduli do not comprise a substantial fraction of the dark 
matter, we require that p$j s be less than the value for baryons, ps/s ~ 10 -10 GeV. If oscillation begins when 
H(T osc ) = m,p, then <f> overdoses the universe for > 10 -26 eV |29j . Since we expect ~ m.3/2, there are 
no viable models with such low scales of SUSY breaking. 

If oscillation begins at earlier times, however, the p$ per comoving volume is suppressed and larger <p masses 
are allowed. One example of how this could happen is if the negative terms proportional to H 2 in the moduli 
potential ((57)) arise only due to couplings with infiaton. If this is the case, the mass-squared of <fi is driven 
positive when the infiaton decays and oscillation begins at the time of reheating, giving 

P4> m l M p\ „ Pb 



< ^. (115) 



For example, if = 0.1 keV and Trh = 10 10 GeV in the scenario where oscillation begins at reheating, the 
present-day moduli density is smaller than the baryon density. 

3.3.6 Mixing due to flat directions 

Combinations of fields with vanishing T— or 2?— terms are also flat and can have large VEVs in the early 
universe. Such directions are typically lifted by higher-dimensional Planck-suppressed operators. For example, 
the potential of a flat direction lifted by a superpotential term of dimension 4 is 

y(^) = (m 2 -a 2 J ff 2 )|0| 2 + -^|0| 6 . (116) 

The interplay between the negative mass-squared coming from the background energy and the <j) 6 lifting term 
gives (j> a VEV for H > m/a, 



M^aH 2 -m 2 y 1/4 



3c 



(117) 



For H 3> m/a, (<fr) ~ H. This is much smaller than the moduli VEV, which is typically Mpi. As a result, 
mxN is also suppressed relative to the moduli-induced mass mixing and the asymmetry transfer is out of 
equilibrium since it was only marginally so for much of the parameter space with mass mixing. The outcome 
is similar to the case for moduli when T\ cp <C T c and ttixn was suppressed. This favors TcV scale or higher 
dark matter. 

There are additional complications to flat direction VEVs giving mass mixing. The direction coupling to 
XN cannot be comprised solely of MSSM fields, since the lowest allowed lifting operator would be dimension-4 
in the superpotential and lead to rapid dark matter decay. The flat directions can also induce mass mixing 
between X and other fields in the theory through renormalizable superpotential couplings, giving X a very 
large mass and lifting it out of the theory. Such models typically require more complicated field content and 
symmetries to give large asymmetry transfer from N — > X . Although possible, the smaller field values make 
models more cumbersome so we focus on the simpler cases in this paper. 

3.4 Mixing induced by background energy 

An alternate, even more generic origin of mass mixing at early times is through couplings to fields dominating 
the energy density of the early universe. We consider the specific case of couplings to kinetic terms of relativistic 
fcrmions i/i in a thermal background. The higher-dimensional operators are 

AC D -4- (#V£>^) (XN + h.c.) . (118) 

In a supersymmetric theory, such terms arise from higher-dimensional corrections to the Kahler potential. If 
%j) are fields dominating the background density, we have 

= ^9*T\ (119) 
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where <?* is the number of relativistic degrees of freedom. We have a mass mixing between scalars 



2 n 2 cg* T 4 tt 2 c 2 
m ^ = ^0M|r 82^^- (120) 

With a higher-dimensional coupling (|118|) . the decay rate is suppressed by (mx/Mp\) and for mx ^ 10 6 
GeV, the dark matter is stable. 

Unlike the case of moduli-induced mixing, direct mixing of fermions by a kinetic term occurs only with 
operators of higher dimension than those inducing scalar mixing. We show below that the asymmetry transfer 
mechanism between scalars is viable only for very high temperatures (around the GUT scale), and this would 
be pushed even higher if the mixing were suppressed by an additional power of T/Mj>\. Therefore, it is only 
the scalar mixing through the thermal background that is relevant. 

We first estimate the expected parameters that give rise to a viable model and subsequently solve the 
density matrix evolution equations numerically. The total thermal interaction rate, averaged between X and 
N, and using the couplings from Section 13751 is approximately 

r.~i(/ + i/)r, (i2i) 

where y and g are Yukawa and gauge couplings given in (1901) . Because the mass mixing is proportional to H, 
it falls as T 2 , and therefore asymmetry transfer is highest when T is large and leptogenesis occurs at a high 
scale. Consequently, the out-of-equilibrium conditions and neutrino masses are consistent with y ~ 1. We 
therefore expect the hierarchy yT ^> mxN to hold. The mass splittings and mixing angle in this limit 
are 

e « y -T, (122) 
6 

sin0 « (123) 

For g ~ 1, it is always true that To ^ e, and so pf?)) reduces to 

2tt 2 T /m X N\ 4 



V-vXt = :pr sin2 cos2 ' 



(^) • d24) 



T 9( ff 4 + y 2 /2) V T 

Integrating over a Hubble time and substituting (|120[) gives 

E^„.„ 7 ^(0. (125, 

For temperatures around Tqut ~ 10 16 GeV, we obtain a fractional X asymmetry (relative to lepton number) 
of around 5 x 10 -4 , which would be the correct asymmetry for a few TeV dark matter candidate. In reality, 
the accumulated asymmetry is larger by a factor of ~ 5 when we account for integrating over multiple Hubble 
times, as we find when we solve the full density matrix evolution equations. 

We solve the system of equations (|55p numerically. We present three benchmark points giving f2rjM ~ 5f2b 
with y = 1 and c = 10: 

• Ti C p = 1.8 x 10 16 GeV, g = 1 and m x = 1 TeV, 

• T lcp = 1.8 x 10 16 GeV, g = 3 and m x = 10 TeV, 

• Ti cp = 5 x 10 16 GeV, g = 10 and m x = 50 TeV. 
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We have assumed that higher dimensional operators mix X with all three generations of N . Asymmetries 
consistent with higher mass dark matter particles are also possible when T\ ep < 10 16 GeV or for g < 1, but 
in these cases, the couplings are insufficiently large to eliminate the symmetric component. Asymmetries 
consistent with lower masses of X are also possible for higher reheat temperatures, but since Trh is already 
borderline given the gravitino problem and the risk of reintroducing the flatness and monopole problems |30j , 
we do not consider any higher reheat temperatures. 

The temperature scales required for scalar mass mixing are higher than we may typically expect from a 
supcrsymmetric theory due to the gravitino problem. The gravitino problem can be alleviated in theories with 
very heavy gravitinos, where the gravitinos decay into axions (31 j or some other non-Standard Model particle. 
Alternatively, one may envision a non-supersymmetric theory with mixed scalars whose scales are stabilized 
by some unknown mechanism. 

Note also that the asymmetry transferred by mass mixing is dominant over transfer from thermal scatterings 
from (|118p . The rate of such processes is 

Tx^nz ~^<H. (126) 

Integrating the rate over a Hubble time gives T 3 /Mp[ ~ 10 -7 even for the highest temperatures we consider, 
T ~ 10 16 GeV and is therefore subdominant to the effect of mass mixing. 

3.5 Phenomenology 

Unlike the two Higgs model in Section [2] the minimum interactions necessary between dark matter and the 
visible sector for the models in this section are Planck-suppressed and there is no built-in mechanism for the 
annihilation of the symmetric component. As a result, the phenomenology is model-dependent. The main 
hope of seeing such a model is through the mechanism for annihilation. In Scctions l3.3l and l3.41 we considered a 
U(l)' gauge interaction under which X and ./V are charged, although other possibilities exist. Direct detection 
bounds strongly constrain the strength of this interaction with quarks. As a result, signals would be small at 
hadron colliders. 

As in Section f2.6l we do not expect asymmetric dark matter models to have large indirect detection signals. 
The mechanisms discussed there, however, could yield a non-zero anti-dark-matter population at late times, 
and so the model would not necessarily be ruled out by a confirmed indirect signal. 

Since our model consists of X mixing with a lepton-number-carrying gauge singlet, as is present in Dirac 
leptogenesis, the determination of the Dirac or Majorana nature of the neutrino masses can give insight into 
the components of our model. However, if neutrinos are discovered to have Majorana masses, this will not 
rule out the models listed above but would favor a more complex model than the minimal Dirac scenario or 
suggest Affleck-Dine mechanism for generating original asymmetries. 

4 Conclusions 

The global symmetries of the universe could have been very different in the early universe than they are today. 
In particular, independent dark matter and baryon asymmetries at late times can emerge from a common 
symmetry at early times and explain the baryon-dark matter near-coincidence. In this paper, we show how 
large field backgrounds in the early universe can induce mass mixings between dark matter and visible matter, 
sometimes leading to distinctive relations between the dark matter and baryon densities that accommodate a 
wide range of dark matter masses. Such mixing can arise in renormalizable Lagrangian terms with field VEVs 
linked to scales already present in the theory, such as the scale of electroweak symmetry breaking. In this 
case, the dark matter mass is naturally either ~ GeV or in the range 300 — 500 GeV. More generically, mixing 
can also arise from Planck-suppressed operators, which we expect to break the global symmetries present 
at low scales. This leads to mixing even in scenarios where there is no renormalizable connection between 
dark and visible matter fields. In this case, the transferred asymmetry is consistent with dark matter masses 
above 1 GeV, with weak scale masses preferred in the case of moduli-induced mixing and leptogenesis at the 
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intermediate scale or in the case of mixing due to background energy when leptogenesis happens at the GUT 
scale. The dark matter masses are typically higher for lower leptogenesis scales, with the only upper bound on 
the dark matter masses being the requirement that the symmetric component be annihilated by perturbativc 
couplings. 

We summarize our results for the various cases we consider in the table below. 
• Two-stage phase transition: 





Relativistic X 


Thermal-dominant X 


Mixing-angle-dominant X 


Fermion DM mass 


O(GcV) 


300-500 GcV 


400-500 GcV 



• Moduli-induced mixing: 





Heavy moduli 


Light stable moduli 


DM mass with 
scalar mixing 


T RH ~ 10 8 GcV: 90 GcV-100 TcV; 
Trh ~ 10 10 GcV: ~ O(GcV) - 
100 TeV 


Transfer highly suppressed 


DM mass with 
fermion mixing 


Trh ~ 10 8 - 10 10 GeV: O(GeV) 


T RH ~ 10 9 GcV: - O(GeV)-100 TeV; 
Trh ~ 10 11 GeV: ~ O(GeV)-100 TeV 


Background energy induced mixing: 



DM mass with scalar mixing T RH ~ 10 ib GeV : 1-100 TeV 



This work expands the range of possible transfer mechanisms in asymmetric dark matter models and 
provides new considerations for dark matter model building. We see that wide range of dark matter models 
and broad range of masses is possible. It is important to explore such possibilities and possible means to 
distinguish among these and other models so we can eventually learn the true nature of dark matter. 
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Appendices 

A Two Higgs model vacua 
A.l Vacua and stability 

The vacua of the potential can be determined by finding the extrema of the potential. Here, we follow the 
analysis of [10| . We perform the analysis at finite temperature and then discuss the particular case of the zero 
temperature vacua. 

Extrcmizing the potential gives 



W=0, T 1 , (127) 





- 2fc 3 (0) 2 


-aiT 2 


fei 


<4l4 


- 2fc 3 (/ l ) 2 


-a 2 T 2 
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At a given temperature, there are four distinct critical points of the potential: 
1. 



(h) = {4>) = 0, (129) 



2. 



W=0, W^i^, (130) 




(0>=O, W = W ^ , , (131) 



4. 



... 8k 2 (4 - 4k 3 fi 2 + (a 2 k 3 - 2a 1 k 2 )T 2 8k^ 2 - 4k 3f i 2 + ( ai k 3 - 2a 2 k 1 )T 2 

{h) -^ 8hk 2 -2kl ' (</,) -\/ 8hk 2 -2kl ' (132) 

For extrema 2 and 3, the vacuum is only well-defined if 

T 2 <T* = ^i. (133) 

Oil 

These temperatures are sufficiently important that we label them critical temperatures and they indicate the 
point at which the fields h and <j) can condense. 

We now examine the stability of the extrema of the potential. 

1. The symmetric vacuum is only stable for T 2 > max(T c 2 1 , T 2 2 ). Below these temperatures, the negative 
mass-squared parameter in the Lagrangian dominates and tends to drive the fields to non-zero values. 
At high temperature, the scalars begin in the symmetric phase and later condense when T passes below 
one of the threshold temperatures. 

2. The vacuum where (h) = and (0) ^ is stable given T < T c2 and k 3 a 2 (T 2 2 - T 2 ) > 2fc 2 ai(T c 2 1 - T 2 ). 
For T c2 > T > T c i, this condition is automatically satisfied. If T < T c \ and T c2 , then the vacuum is only 
stable down to the temperature 

2 _ 4(2^4-^4) 
k 2 a.\ — k 3 a 2 

at which point this cxtremum becomes a saddle point. 

3. The vacuum where (0) = and (h) ^ is exactly the same as the one in vacuum #2, except with indices 
1 <-> 2. 

4. Finally, the conditions for a stable vacuum at both h and <p non-zero are 

2fc 2 a 1 (T 2 1 - T 2 ) - k 3 a 2 (T 2 2 - T 2 ) > 0, (135) 
2fcia 2 (T 2 2 - T 2 ) - ha^T 2 , - T 2 ) > 0, (136) 
4fcifc 2 - fcf > 0. (137) 

It is not immediately apparent whether this is stable or not, but we demonstrate below out that for all 
parameters of interest to us (namely those that give a two stage phase transition), the conditions are be 
violated and this is be a saddle point. 
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Because externum # 4 is a saddle point, there is only one non-zero VEV in any vacuum. This gives us two 
options: the first is that the field that initially condenses is also the field with a VEV at zero temperature (we 
want this to be h) and hence there is only a one stage phase transition. The other possibility is that the initial 
field to condense (<p) has a higher zero-temperature energy when condensed than the vacuum with non-zero 
h, and so at late times there is a phase transition where the <fi VEV shuts off and the h VEV is turned on. 
Because there is a classical energy barrier between these two minima of the potential, this generally leads to 
a first-order phase transition (see Appendix |B|) . 



A. 2 Parameter constraints 

In this section, we determine the parameters that lead to a viable two stage phase transition. We begin by 
matching the parameters onto the known zero temperature physics. In particular, we would like {<p) = and 
(h) = v at zero temperature, where v = 246 GeV. To ensure that h condenses instead of <j), we require the 
vacuum energy in the state where h is condensed to have a lower energy than the state where <f> is condensed. 
The vacuum energy in the (cf>) — vacuum at zero temperature is 

Vi(T = 0) = -£ (138) 
fei 

and similarly for the (h) = vacuum energy Vi. For the former to be lower, we require |Vi(T = 0)| > |V2(T = 
0)| or, alternately, 

ti > h 

Imposing this constraint, and therefore being in the (</>) = vacuum at zero temperature, we can determine 
the physical Higgs mass from the quadratic term, 

ml = 8fil (140) 

The zero temperature Higgs mass has not yet been determined, but electrowcak fits prefer a value as close as 
possible to the LEP bound of 115 GeV, so we choose tuh = 120 GeV. We can then determine the parameters 
fii and fei using 

/tl = r, 42.4 GeV, (141) 
V8 

h = l(-) 2 w 0.0297. (142) 



4 > 7 1 - ( 139 ) 



V V 

Meanwhile, <f> gets a mass due to the 4> 2 h 2 term, and its value is 



This gives us our second constraint, 



m 2 = 2^-4/4 (143) 



k 3 > ^4h, (144) 
Mi 



which requires that the mixed quartic term be large enough to compensate for the negative <p mass-squared 
term. If ()144[) is not satisfied, then <f> is also driven to a non-zero value (vacuum #4 above). 

In order to have a two stage phase transition, we would like the system to first condense in the <f> direction. 
This imposes the requirement that T c \ < T C 2, which can be expressed as 



4 < ^; d45) 
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combining this with (|139p gives 




(146) 



Thus, we see that the mass terms are constrained to be roughly close in magnitude to one another in order to 
have a two stage phase transition. 



B The two Higgs phase transition and bubble nucleation 
B.l Tunneling and bubble nucleation 



There exists a potential barrier between the true and false vacua in our two Higgs model. The vacua become 
degenerate at the temperature 

4 (y/kin\ ~ Vh^l) 

l -d ~ 



Ta = 



(147) 



/feCKl — Vfcl a 2 

Below this temperature, the system can transition from the false vacuum to the true vacuum. The phase 
transition between the two vacua, should it occur, satisfies the following conditions: 

1. It occurs at the very latest before big bang nucleosynthesis (T 2> 10 MeV) to satisfy observational 
constraints. 

2. The 4> VEV must turn off sufficiently rapidly that the asymmetry relation from mass mixing is not 
disturbed. This means either that the fields directly tunnel into the = vacuum, or in the case where 
the system tunnels to a non-zero <j> vev, that <j) is rapidly damped to zero (on a time shorter than the 
rate of asymmetry transfer) following tunneling. 

3. The expansion of the bubble wall should be faster than the asymmetry transfer rate from thermal 
scattering. 

To begin, we calculate the temperature at which bubble nucleation occurs. At large temperature, when 
the bubble wall r w 3> /3, we can take the bubble solution to be constant in the period Euclidean time, r. The 
bounce solution is then given by the 0(3) symmetric function minimizing the Euclidean action 



S{T,ct>,h) = 



4tt 



dr r 



dr 



V(<f>,h,T) 



(148) 



The semi-classical probability P of nucleating a bubble inside of a casual volume during the temperature 
interval dT is [32] 



(149) 



where a; is a dimcnsionlcss coefficient given by functional determinants from the path integral and £ sa 1/34 is a 
numerical factor coming from the Hubble constant. The result is dominated by the action S in the exponential, 
making the precise value of uj irrelevant. Taking to ~ 1, T ~ 30 GeV and AT ~ 10 GeV (for reasons that 
will be shown later), we find that S « 140 gives a probability of nucleation equal to one. We take this as our 
benchmark for determining when nucleation occurs. 

It is difficult to find the bounce solution to (|148[) numerically since there are two fields involved. Because 
we are simply interested in an estimate of the action, we instead make an ansatz that eliminates one of the 
fields. We take h = /(</>), with / to be determined from V{4>, h), giving 



S(T, 



T 1 



V(4>,T) 



(150) 
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The resulting equation of motion is 



d 2 (f> 2 d(j) _ 1 
'ch J+ r~ch : ~ l + /'(0) 2 



dV 



(151) 



The initial condition 0(0) is chosen to satisfy the bounce boundary conditions 4>'(0) = and 4>{oo) = 0f a ise- 
This equation can be solved numerically using the "overshoot-undershoot" method [33] . 

The ansatz h = f((f>) is chosen as the curve in (f> — h space such that the smallest possible potential barrier 
is crossed while tunneling. We determine this curve by slicing the potential along lines of constant <j> and 
determining the value of h that minimizes V. While we cannot guarantee that this is the minimum action, 
it is plausible that it is among the curves that gives the smallest action, and it gives an upper bound on the 
action. If we find rapid bubble nucleation from our ansatz, then bubble nucleation does occur, even if the 
exact instanton configuration is slightly different from the one determined here. 

Bubble nucleation does not occur for all choices of parameters; for a particular choice of k^, hi must be 
chosen sufficiently large to reduce the barrier between vacua. Numerical solutions are presented in Section [2 .41 

To address point #2 above, we find that tunneling does not, in fact, take <f> to zero but instead to some 
non-zero value. then undergoes damped oscillations about its minimum, = 0. The damping due to the 
Hubble constant, H ~ 10 -13 GeV, is negligible. The dominant source of damping is decays of — >• X + L 
and the conjugate process. These decays do not change the asymmetry in the X or L sectors, since it is 
CP-conserving. We require that r^-s.^ to avoid eliminating the asymmetry in X by the residual mass 

mixing due to non-zero 0. We can estimate by 




(152) 

The asymmetry washout term is [71 with Tq given by the larger of 

4 2 

T x-^^§T 3 , (153) 



ni 



r L ~ j|r. (154) 

The first interaction comes from XX — > LL through an intermediate 0, while the second is L interacting with 
W bosons, whose mass is given dominantly by the h VEV. For a benchmark point with yx = 1, mx = 400 
GeV, m<f, = 450 GeV and Tn = 60 GeV (as is a typical scale for bubble nucleation), we get = 16 GeV while 
ro ~ r^, = 7.6 x 10~ 2 GeV. Thus, the VEV of <j> safely damps away before the asymmetry relation from mass 
mixing is altered. 

Finally, the bubble wall must sweep over a particle without it undergoing any thermal interactions. There 
is approximately one particle in a region with length 1/T. The bubble wall takes a time t w = (2Sr + l/T)/v w 
to pass over the cube, where 8r is the bubble wall thickness and v w is the velocity. This time must satisfy 
t w <C 1/Tl-^x- We find that the wall thickness is always < 1 GeV -1 . Using the rates of interaction in the false 
vacuum phase and assuming a small mixing angle, we determine that v w > 10 -3 for the bubble wall passage 
to be out of equilibrium. This is easily satisfied for variants of the Standard Model with a first-order phase 
transition (the drag from top quarks in that case is similar to the drag for this bubble wall) [32]. Therefore, 
the transition is strongly first-order and our instantaneous approximation for the state projection is valid. 

For GeV-scale dark matter, where the mixing angle is large, there is no thermal suppression of any of the 
modes (n X i = til 1 ) and neither the number densities nor the asymmetries change during the phase transition. 

B.2 Back-reaction on phase transition 

So far, we have neglected the back-reaction of the mass terms on the evolution of the bubble expansion at the 
first-order phase transition. In the limit of T t <C M, the dark matter energy density present after the phase 
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transition is entirely given by the repopulation of the heavy states due to the mass mixing. We must also 
concern ourselves with the "drag" associated with the top quarks becoming massive at the phase transition. 

The energy density of the light mass eigenstates that become heavy X fields following the phase transition 
is given by 

2 

px = Tax sin 2 6riL< — — ^ sin 2 OmxT™, (155) 

while the mass energy density of the top quarks following the phase transition is 

12 ( v h (T N ) \ 3 
Pt = mtnt = ^{2l6Gcv) mtT v (156) 

The contribution from top quarks dominates over the contribution from X even for the largest viable masses 
for our model mx ~ 500 GeV, since there are six top quark degrees of freedom, giving an "effective" top quark 
mass of 6m t ~ TeV. 

At high temperature, the quarks have a large kinetic energy and this can be used to offset the vacuum 
energy expended to impart mass to the quarks. The energy density from the kinetic energy of the quarks 
(using (p) ps 3T) is 

Pk = (157) 

7T 

and so for T t > 50 GeV, we find that the mass of the top quarks can be accounted for entirely from the kinetic 
energy and so this does not affect the progress of the phase transition. 

In the region 30 GeV < T t < 50 GeV, the situation is slightly more subtle. The energy density that goes 
into creating masses for the top quarks is comparable to the energy difference between the vacua immediately 
after T t . We must therefore not use T t as the tunneling temperature and use instead the temperature at 
which the energy difference between the vacua is equal to the energy density from top quark masses. For a 
benchmark point with yx — 0.4, = 0.5, £3 = 0.7 and p2 = 80 GeV, we find that Tn = 34 GeV, whereas the 
temperature where the energy difference between vacua is pt is T = 30 GeV, corresponding to an approximately 
10% difference. Similar results are found for nearby benchmark points. 

Bubble nucleation rarely occurs below 30 GeV, so we do not consider this range further. 

Fortunately, this does not significantly alter our results, since in the mass mixing-dominated limit, the final 
X density is fixed only by tan 2 8, which has only a weak dependence on T. In fact, this gives us slightly more 
flexibility as the temperature for bubble nucleation can be about 10% lower than before, giving more thermal 
suppression and a broader range of parameter space where the X density is dominated by mass mixing. 



C Review of Dirac leptogenesis 

We review a mechanism for realistic Dirac leptogenesis [19) . As the universe cools below the temperature 
T = -My,, ip and ip decay in two ways: 

i> N^Hl (158) 

i, -+ L X - (159) 

$^NH U , (160) 

rp -> LV- (161) 

The total decay rates of ip and ip are equal, but they can decay with different rates into N and L. This can give 
rise to an asymmetry of L over . Because total lepton number is conserved by these interactions (assuming 
ijj carries L = 1), this asymmetry is equal and opposite to the asymmetry of TV over iV T . Thus, there is an 
iV asymmetry that is directly correlated with the Standard Model lepton asymmetry (and hence related to 
the baryon asymmetry by sphaleron processes) . N and L ultimately come into equilibrium due to oscillations 
once H u and x have VEVs (see [T2"]). 
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Eventually, TV and L come into thermal equilibrium and the lepton asymmetry is destroyed at the time of 
cquilibriation. This must happen after the time of BBN, since strong constraints exist on the Standard Model 
lepton asymmetry in that era. Fortunately, this equilibriation rate is suppressed by neutrino masses and is 
not important until T ~ eV. We can see that this is true below by integrating out ip, giving an effective 
superpotential 

W e g D MEMK NHuLx . (162) 

At low temperatures, \ and H u have weak scale VEVs (which come about due to interactions with the SUSY 
breaking sector), giving neutrino masses 

m v ~ — < eV, (163) 

meaning that ul Vn ~ 0.1 — 1 for ~ 10 15 GeV and weak scale (x), while it is i/l Vn ~ 10 -5 for M$ ~ 10 10 
GeV. 
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